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ABSTRACT 

The outbursts of low mass X-ray binaries are prolonged relative to those of dwarf 
nova cataclysmic variables as a consequence of X-ray irradiation of the disc. We show 
that the time-scale of the decay light curve and its luminosity at a characteristic time 
are linked to the radius of the accretion disc. Hence a good X-ray light curve permits 
two independent estimates of the disc radius. In the case of the milli-second pulsars 
SAX J 1808.4-3658 and XTE J 0929-314 the agreement between these estimates is very 
strong. Our analysis allows new determinations of distances and accretion disc radii. 
Our analysis will allow determination of accretion disc radii for sources in external 
galaxies, and hence constrain system parameters where other observational techniques 
are not possible. We also use the X-ray light curves to estimate the mass transfer rate. 
The broken exponential decay observed in the 2002 outburst of SAX J 1808.4-3658 
may be caused by the changing self-shadowing of the disc. 

Key words: accretion, accretion discs, binaries: close, x-rays: binaries, stars: indi- 
vidual: XTE J 1808.4-3658, stars: individual: 4U 1705-44 



1 INTRODUCTION 

King & Ritter (1998) (henceforth KR) and Shahbaz, Char- 
les & King (1998) (henceforth SCK) examined the X-ray 
light curves of transient low mass X-ray binaries (LMXBs) 
in decay, characterising them as simple exponential or linear 
decays. In some cases, however, the decays show a distinct 
knee. As shown in Fig. 2, this is not a secondary maximum 
in the sense that the X-ray count rate declines throughout. 

In X-ray bright disc-accreting systems, including tran- 
sient LMXBs during outburst, the temperature of the ac- 
cretion disc is dominated by X-ray heating from the inner 
accretion regions across most of the disc. For a disc in which 
the scale height can be described as a function of radius by 
H = HoR n for constant Ho and n, the temperature, T, is 
given by de Jong, van Paradijs & Augusteijn (1996) as 



1-AH 



(n-l)Lx, 



(1) 



AiraR 2 R 

where A is the disc albedo, Lx is the central X-ray luminos- 
ity heating the disc, assumed to be emanating radially from 
a small spherical source. Only if n > 1 is the whole surface of 
the accretion disc illuminated; if n — 1 (1) predicts a surface 
temperature of zero which is invalid because of local viscous 
heat production. Consistent with KR we adopt a value of 
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n — 9/7, although see section 5. (1) can be rewritten to give 
the maximum radius Rh which is heated by X-rays to the 
temperature needed to remain in the hot, viscous state. 
Thus 
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jH (n-l)L x . 



We abbreviate this equation as 
Rt n = $L X . 



(2) 



(3) 



We note that especially in the case of a black hole the X-ray 
source will be the inner accretion disc, and as such disc-like 
rather than point-like. This introduces an factor of H/R into 
the r.h.s. of (2), which is equivalent to modified values of n 
and $ in (3); specifically, the index of 3 — n in the l.h.s. 
becomes 4 — 2n, and $ gains a factor of Ho- However, n 
and Ho may differ systematically between neutron star and 
black hole systems, in which case the differences in (3) will 
be less simple. 

KR showed that X-ray heating during the decay from 
outburst causes the light curves of transient LMXBs to 
exhibit either exponential or linear declines depending on 
whether or not the luminosity is sufficient to keep the outer 
disc edge hot. They note that exponential decays must re- 
vert to the linear mode when the X-ray flux has decreased 
sufficiently, but do not analyse this in detail. Nor do they 
consider the effects of mass transfer from the donor star, 
— M2, during the outburst. Whilst in many systems — M2 
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is negligible compared to the mass accretion rate onto the 
compact object, M c , during outburst, this is not necessarily 
the case. In this paper we consider these issues. Section 2 
contains our theoretical predictions for the decay light curve. 
In section 3 we identify some suitable light curves and de- 
scribe their parametrization. We compare the results from 
the fitting in section 4. 



2 ANALYSIS 

2.1 Exponential decay 

The mass of the accretion disc during the exponential decay 
is given by equation 3 of KR; 



M di! 



(4) 



where vkr is some measure of viscosity, taken by KR to 
be the viscosity near the outer disc edge with a value of 
fKR ~ 10 11 m 2 s _1 . In the case of significant mass transfer 
into the disc from the donor star the central accretion rate 
is given by 



M c = -M 2 - M disc . 

Using this to replace M c in (4) we obtain 
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Therefore the mass of the hot disc can be written 

. .. ..... . #disc(-^) 

M disc = M a exp [ — ) + ^ '-, (7) 
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where M a is the constant of integration, and the X-ray lu- 
minosity is proportional to 



M c = -M 2 H -J exp 

^disc 
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At some time, tt, the temperature of the outer disc edge, 
while still dominated by X-ray heating, will be only just suf- 
ficient to remain in the hot, viscous state, i.e. Rh = -Rdisc- 
We denote the corresponding X-ray luminosity as L t , with 
the central accretion rate M t . The X-ray luminosity at ear- 
lier times is 



L x = (L t -L 2 )exp|- ^#^l )+L 2 , 



r>2 
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(9) 



where L 2 = r) ( — M 2 j c 2 , and r\ ~ 0.1 is the efficiency with 
which rest mass is liberated as X-rays. 



2.2 Linear decay 

When Rh < -Rdisc the outer part of the disc is no longer 
kept in the hot, viscous state by central irradiation. The 
radius of the hot part of the disc will decrease with the 
decreasing X-ray luminosity, corresponding to the late linear 
decline of KR. This, too, may be modified by the donor 
star mass loss rate, depending on the mass accretion rate 
through the now cold outer disc. Designating /i c (R) as the 
rate at which mass is transported inwards through the cold 
disc, we examine the two extreme behaviours of fi c . For the 
minimally efficient cold disc /i c (R) = for all R, and none 



of — M 2 reaches the hot disc when Rh < -Rdisc. At the other 
extreme /i c (R) = — M 2 and the cold disc does not increase 
in surface density with time. The occurrence of outbursts 
suggests that the latter case is false; with stable cold-state 
mass transfer, outbursts would never occur. 

The mass in the hot disc changes due to three factors: 
accretion from its inner edge, the loss of material into the 
encroaching cold disc, and mass transfer from the cold disc. 
Hence, 

Mhot = -Mc + Vc(Rh) + 27rRfcS(R h ) R h . (10) 
Using the approximate substitution from (3) that 
Rl = <S> V c 2 M c , (11) 
(10) becomes 

Mhot = -M c + ti c {R h ) + 7r<I> J? c 2 S(R h ) M c . (12) 
Additionally, equation 2 of KR allows Mhot to be written as 
M c Rl 



M hot = „ 

Equating (12) to d/dt (13), 
1 d 



(13) 



M c Ri = -M c + n c (R h ) + 7r$77c 2 S(-R h ) M c , (14) 



3i/ K r dt 
so using equation 2 of KR 

^ = 3i/ K R ( Mc(Rfe) _ 1 
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We now substitute the two cases of —\i c (Rh) outlined above 
When mass transfer into the hot disc is negligible, 

3^kr 



M, 



$7/C 2 ' 



(16) 



and the decay is the linear decline predicted by KR as ex- 
pected. When fic(Rh) = — M 2 , 



3z/kr (-Mz - Mc) 



M c 



(17) 



While M c S> — M 2 this reduces to the same form. As M c 
decreases it will approach the limit of M c = — M 2 , at which 
point M c = 0. The exact form of the decay will depend on 
the unreliably known form of ^i c (R), and during the late lin- 
ear decline viscous heating near Rh becomes non-negligible. 
We do not predict the form of the late decay other than to 
note that if there is some radius Rn m satisfying 



/i(Rlim) f]C 



(18) 



then the l.h.s. of (18) will give the asymptotic limit of the 
decay and the accretion luminosity at the beginning of qui- 
escence. 

As a final point on the linear decay we note that the 
viscous instability may cause a range of annuli within the 
cold disc to enter the hot state and rapidly transfer mass in- 
ward. Some of this may enter the inner hot disc, or else may 
increase the surface density of the subsequent cold disc such 
that the inner hot disc receives more material at its outer 
edge. Consequently, small rebrightenings or other artefacts 
potentially complicate the linear decay. 
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2.3 Continuous derivative 

To examine the observed X-ray light curves of transient 
LMXBs it is useful to examine the constraints on the tran- 
sition. The gradient of the exponential decay at time t t is 

3v K r 



Lx = ~ 
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(L t - r, (-M 2 ) c) 



(19) 



whilst the gradient of the linear decline is given by (16) 
3^r 



(20) 



From (3) we adopt 



(21) 



where it is assumed that at the outer disc edge H ~ 0.2i?dh 
We therefore obtain 



Lx(U) 
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(22) 



In the case that — 77M2C is small relative to L x the gra- 
dient of the exponential decay at time t t is equal to that 
of the subsequent linear decline; the first derivative of the 
X-ray light curve is expected to vary smoothly throughout 
the decay. If however — 7/M2C 2 cannot be neglected then the 
linear decline is steeper than the exponential decay at t t and 
a discontinuity in gradient is expected. Whether this is de- 
tectable depends on the value of — M2 and the quality of the 
X-ray data. 

Some X-ray light curves, including most of those dis- 
cussed in this paper, include a knee feature in the decay 
from outburst. Where this can be explained by an exponen- 
tial to linear transition with non-negligible —r\M2<? we label 
it as a 'brink'. In general the decay light curve may contain 
multiple knees due to different effects, but we do not expect 
more than one brink unless the source rebrightens between 
them such that each occurs at approximately the same lumi- 
nosity, corresponding to a disc not varying greatly in radius. 



3 APPLICATION TO OBSERVED X-RAY 
TRANSIENTS 

X-ray light curves from the All-Sky Monitor (ASM) and 
Proportional Counter Array (PCA) of the RXTE space tele- 
scope were examined to find examples of decay light curves 
which fit the template predicted above. We use our analy- 
sis to determine distances and accretion disc radii. System 
parameters from the literature for the objects analysed are 
collected in Table 1 for comparison with quantities we de- 
rive. Since we will be calculating the radius of the accretion 
disc we calculate from the system parameters the orbital 
separation a. 

The count rates N as a function of time t were fitted 

with 

N = (N t — N e ) exp {- t -^r 1 ) + Ne (23) 

for the exponential phase, where N e is the limit of the ex- 
ponential decay, Nt is the count rate at the brink, and r e is 
the time-scale of the decay. The linear decline is given by 

t-t t - 



N = N t 



(24) 



where ti is the time after the brink at which the count rate 
would become zero if the linear decline continued. Using 
the assumed distance to each source with a hydrogen col- 
umn density of 10 22 cm~ 2 it is possible to convert the count 
rates Nt and N e into absolute X-ray luminosities Lt and L e , 
the X-ray luminosity at the brink and the limit of the expo- 
nential decay (corresponding to —M2) respectively. For the 
two systems for which we have simultaneous ASM and PCA 
light curves, XTE J 0929-314 and GRO J 1744-28, the ratio 
of count rates is approximately 31. We assume this conver- 
sion to be valid for all sources for comparison purposes, and 
present luminosities in the 1.5-12 keV range corresponding 
to the ASM spectral range using an assumed power law spec- 
trum with photon index of F — 2. The luminosities thus de- 
termined are given in Table 2. We note that if the spectrum 
differs between systems our deduced Lt is invalid, though 
our deduced exponential timescale will be unaffected. 

Using (3) and (9) it is possible to calculate the disc 
radius independently from the measured values of Nt and 
r e ; the result based on r e is of particular importance as it 
does not depend on the conversion of count rate to flux but 
only on a directly observable timescale. It is also possible 
to determine — M2 from N e . Results are given in Table 3; 
for comparison we have also listed the circularization radius, 
i? c i rc , and the distance between the compact object and the 
inner Lagrange point, b\, for each source based on the values 
in Table 1. Discussion of the choice of values for <E> and vkr 
is given in section 4. 



3.1 4U 1543-475 

Fig. 1 shows that the exponential-then-linear fit to the decay 
from the 2002 outburst of 4U 1543-475 is significantly better 
than the least-squares regression line. Because no disconti- 
nuity is seen in the gradient, the fitting routine used was 
based on the continuous derivative model (see section 2.3). 
Previous outbursts were reported in 1971, 1983 and 1992. 
The source therefore spends only a small fraction of its time 
in outburst, suggesting that — M2 is small compared with 
the central accretion rate during outburst, consistent with 
the source's failure to display a broken gradient. 



3.2 SAX J 1808.4-3658 

SAX J 1808.4-3658 was first detected in outburst in 1996 and 
has subsequently shown outbursts in 1998, 2000, 2002, and 
2005. Hence it is expected that the ratio of — M2 to the peak 
central accretion rate is much higher than in 4U 1543-475. 
This is consistent with the appearance of a broken decay 
in the two well observed outbursts, shown in Figs. 2 and 3. 
Each light curve is well fitted by a decay that is exponential 
towards a positive limit followed by a sharp gradient dis- 
continuity (brink) leading to a linear decline. The fit to the 
1998 outburst includes all points that can be confirmed to 
occur during the decay; the observation at day 16 of Fig. 2 
may belong to a period of constant luminosity or even an 
unsteady rise to outburst, while the observations around day 
33 of Fig. 2 may involve a rebrightening. Even if this is not 
the case the late linear decline is expected to become in- 
creasingly complex, as we noted in Section 2.2. Comparison 
with Fig. 3 shows that the brink occurs at approximately 
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Modified Julian Day - 52 400 

Figure 1. Top panel: The fitted ASM light curve of 4U 1543- 
475. The solid line shows the cxponcntial-thcn-lincar fit, the 
dashed line indicates the best linear fit. The data shown in grey 
were not included in the fit, lying before and after the period 
of interest. Middle panel: Residuals to the linear fit. Bottom 
panel: Residuals to the exponential/linear fit. 
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Modified Julian Day - 50 900 

Figure 2. The fitted PCA light curve of the 1998 outburst of 
XTE J 1808.4-3658. The dashed line indicates the fit of SCK. 



the same count rate in both cases. This is expected if 7?disc 
remains constant (i.e. Lt is the same). 

The dashed line in Fig. 2 indicates the fit proposed by 
SCK, in which they identified the surplus around the brink 
as a secondary maximum. We disagree and view their fit 
as inadequate in having a negative limit to the exponential 
decay. It is possible that features in other decay light curves 
identified as secondary maxima but having monotonically 
decreasing X-ray flux have the same explanation. 
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Figure 3. The fitted PCA light curve of the 2002 outburst of 
XTE J 1808.4-3658 



A second issue in Fig. 3 is the broken nature of the light 
curve prior to the brink. Whilst we do not have a detailed 
explanation for this, we speculate that it may be associ- 
ated with some regions of the disc being self-shadowed and 
initially remaining in the cool state. 3D hydrodynamic sim- 
ulations (e.g. Foulkes, Haswell & Murray 2006) show self- 
shadowing is likely, due to spiral density waves and irra- 
diated disc warping (c.f. Pringle 1986). Consequently the 
effective area of the hot disc, the quantity probably corre- 
sponding to our measurement of R\ iBC , increases during the 
decay in an apparently stepwise manner, introducing new 
material into the hot disc and therefore resetting the cen- 
tral luminosity to a higher value. This view is supported by 
the fact that fitting the region around day 65 of Fig. 3 in- 
dicates a shorter exponential time-scale than from day 70 
to the brink, corresponding to a smaller outer radius of the 
hot region. By day 70 we assume that the entire disc is in 
the hot state. The growing irradiated area could lead to an 
actual increase in L x , hence it is an explanation for a type of 
secondary maximum. It is possible that the light curve will 
monotonically decline despite the growing irradiated area, 
hence producing knees in the light curve. The light curve of 
SAX J 1808.4-3658 around day 67 of Fig. 3 is suggestive of 
this behaviour. 



3.3 XTE J 1751-305 

XTE J 1751-305 was first detected during its 2002 outburst, 
shown in Fig. 4, and a second outburst was observed in 
2005. If this repetition interval is normal for the system, 
it indicates a ratio of — M2 to peak central accretion rate 
more similar to that of SAX J 1808.4-3658 than to 4U 1543- 
475. Consistent with this, the 2002 decay light curve has a 
shape similar to those of SAX J 1808.4-3658 in that a smooth 
exponential decay is broken by occasional rebrightenings. 
Therefore only the last portion, appearing to correspond to 
a single smooth decay, is fitted. The linear portion of this 
decay terminates at a constant level. In accordance with (17) 
we interpret this as the rate at which mass can be transferred 
through the mostly cold disc; this rate will continue during 
quiescence whilst the disc mass increases. 
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Modified Julian Day - 52 370 

Figure 4. The PCA light curve of the 2002 April outburst of 
XTE J 1751-305, fitted after MJD 52 373. 
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Figure 

ASM lie 



5. The fitted PCA light curve of XTE J 0929-314. The 
ht curve, suitably scaled, is overplotted in grey. 



3.4 XTE J 0929-314 

A third observed example of the brink is found in the only 
observed (2002) outburst of XTE J 0929-314. Here the over- 
all shape before the brink is exponential, but a maximum at 
about day 10 of Fig. 5 and a minimum at about day 19 are 
anomalous. Possibly the outer edge of the disc remained cold 
until late in the exponential decay and several secondary 
maxima of the type already discussed have occurred. 



3.5 XTE J 1807-294 

XTE J 1807-294 is, like SAX J 1808.4-3658, XTE J 1751-305 
and XTE J 0929-314, one of seven known milli-second pulsar 
(MSP) systems. MSPs have short orbital periods and hence 
small and therefore relatively simple discs. We therefore ex- 
pect the X-ray light curve of XTE J 1807-294 to also display 
the brinked decay seen in the first three. Fig. 6 does not ap- 
pear to lend itself to this type of fitting, and we have been 
unable to obtain reasonable fitting parameters. One possible 



Modified Julian Day - 52 700 
Figure 6. The PCA light curve of XTE J 1807-294 

explanation is that the linear decay had already begun at the 
start of the observed light curve; in this case the brink count 
rate must be around fifty per cent greater than that seen in 
SAX J 1808.4-3658. The orbital separation of XTE J 1807- 
294 (c.f. Table 1) is half that of SAX J 1808.4-3658, so the 
disc must be correspondingly smaller. A high value of N t can 
only be accounted for if XTE J 1807-294 is at a distance of 
less than about 1 kpc so the luminosity remains low while the 
countrate is high. Conversely, if the light curve corresponds 
to the exponential portion with several rebrightenings, the 
maximum value of N t may be 15, with the light curve after 
day 80 tentatively assumed to be linear. If the disc is half 
the radius of that of SAX J 1808.4-3658, in proportion to 
the orbital separation, (12) indicates that the brink lumi- 
nosity will be one quarter that of SAX J 1808.4-3658. The 
PCA count rate at the proposed brink is one quarter that 
of SAX J 1808.4-3658, indicating a similar source distance 
of D ~ 2 kpc. We deem the latter explanation more likely, 
and suspect that the exponential decay of XTE J 1807-294 
is complicated by self-shadowing of the disc. 

The more extreme mass ratio of XTE J 1807-294 than 
SAX J 1808.4-3658, indicated in Table 1, suggests that the 
accretion disc radius in XTE J 1807-294 is a larger fraction of 
the orbital separation, indicating a larger brink luminosity 
than assumed above, increasing the distance estimates in 
both cases. A lower value of Nt is possible in the second 
case, which would also indicate a greater source distance. 

3.6 GRO J 1744-28 

The earliest part of the RXTE light curve of GRO J 1744- 
28, shown in Fig. 7, shows the system already in outburst; 
the first PCA data are consistent with the outburst decay. 
The light curve has a brink just after day 170 in Fig. 7. 
GRO J 1744-28 entered a relatively quiescent period in 1996 
April (~day 200 of Fig. 7), punctuated by a mini-outburst 
shown inset in Fig. 7 following day 260 which reached less 
than half the brink luminosity. A second outburst, smaller 
than the 1996 outburst, followed during the first four months 
of 1997. Since then the system has been quiescent. Fig. 8 
shows the decay from the 1997 outburst, in which a brink is 
clearly present. Our confidence in the identification of these 
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Modified Julian Day - 50 000 

Figure 7. GRO J 1744-28 light curve for the decay of 1996. Points 
shown in black are the PCA data; circles indicate those data not 
fitted. ASM data is shown in grey to mark the behaviour during a 
gap in the PCA coverage. The solid line indicates the best-fitting 
decay model; the dotted lines indicate the error bars on the model. 
The dashed line marks the best-fitting date of the brink. Inset: 
The ASM light curve for the late decay and subsequent mini- 
outburst. 
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Figure 9. The fitted PCA light curve of GX 339-4. Circles mark 
data points where the y-scale has been stretched by a factor of 
ten. 
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Figure 10. The PCA light curve of XTE J 1550-564. Not enough 
data cover the decay to be able to be able to perform useful fitting, 
but the presence of a knee in the decay is indicated by the four 
points through which lines have been drawn. 



Figure 8. GRO J 1744-28 light curve for the decay of 1997. Solid 
and dotted lines mark the best fit and uncertainty, as in Fig. 7. 



two knees as brinks is increased by their similarity in lumi- 
nosity. 



3.8 XTE J 1550-564 

Although not enough data exist to be able to establish the 
parameters of the brinked decay, a knee occurs in the PCA 
light curve, which has been marked in Fig. 10. Interpreting 
this as the brink it is therefore possible to calculate i?disc 
based on N t but no other system parameters. 



3.7 GX 339-4 

As shown in Fig. 9, GX 339-4 follows an exponential de- 
cay to less than 10 per cent of its outburst maximum, a 
much smaller fraction than is exhibited by the other sys- 
tems considered so far. The relatively few existing points 
fit the brinked decay model well. The late decay departs 
smoothly from the linear decline, in line with expectations 
(c.f. section 2.2 eqn. 15 and following comments). 



3.9 GRO J 1655-40 

The ASM light curve of GRO J 1655-40 shows one outburst, 
starting in 1996 April and persisting until 1997 August. This 
outburst is broadly divided into two broad peaks, of which 
the decay of the second is shown. Although the system had 
been in outburst for a considerable time, Hynes et al. (1998) 
note that no stable hot disc state exists for luminosities 
less than the Eddington luminosity; it is therefore unclear 
whether material from the outer disc edge participated in 
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Modified Julian Day - 50 600 

Figure 11. The ASM light curve of GRO J 1655-40. Although the 
calculated error on the brink luminosity is less than 1 ASM count 
s -1 , we note by examination that the brink in the light curve 
may be uncertain by 5 ASM counts s -1 if some of the scatter is 
intrinsic variability. 



the outburst. As shown in Table 3, the accretion disc radius 
calculated from the brink luminosity and the exponential 
timescale match each other, but are considerably less than 
the system's circularization radius. This may indicate that 
the mechanism examined here is applicable, but acts on a 
small inner disc region rather than on the disc as a whole. 
Possibly this second peak is caused by a high surface density 
remnant left by the preceding outburst activity. 

GRO J 1655-40 returned to outburst in 2005; if the av- 
erage flux between the beginning of the two successive out- 
bursts is taken as an indicator of —M2, this value is found 
to lie below the brink by a factor of three. The mass transfer 
rate from the secondary star may be enhanced during out- 
burst, but it is likely that the mass flow onto the (constant 
radius) outer edge of the hot disc during the exponential 
phase is supplied from additional disc structures rather than 
from the secondary directly. The complexity of the disc in 
this system does not necessarily invalidate the application 
of this model during at least some phases of the decay. 



3.10 4U 1705-44 

No PCA light curves for 4U 1705-44 are available, so the best 
evidence for our model being applicable to this system comes 
from the ASM light curve. This shows numerous outbursts, 
the decays of which are shown in Fig. 12. A range of smaller 
outbursts also occur which are not shown in Fig. 12. Arrows 
mark features in the light curve that we have interpreted as 
the brink of each decay; using these as part of the initial 
conditions for finding a best fit we obtained the best-fitting 
decays overplotted. Clearly there is substantial variation in 
the height of the brink, which indicates variation in the ra- 
dius of the disc. In order to explain this we used the ASM 
light curve to assess the fluence of each outburst. A greater 
fiuence suggests a greater initial disc mass, which leads us 
to expect positive correlation between the fluence and cal- 
culated disc radius. Relative fluences are given in Table 3, 







a 


b 




d 








" '''Vjl 



















2 keV 




f 


"i 1 ' 


9 


l h 


7 20 




- ^ W 


,,, v. 


„'V| 

I 1 J 





1200 1300 1350 




1750 1800 



2050 2100 2150 2200 225': 



2500 2550 2800 2850 2900 3150 

Modified Julian Day - 50 000 

Figure 12. The ASM light curve for all decays from major out- 
burst seen in 4U 1705-44 by the instrument to date. Arrows marks 
points in the light curve where a knee between an exponential de- 
cay and a linear one may be present, assigned by eye and used 
as initial conditions for the fitting. Horizontal dotted lines indi- 
cate the limit of the exponential decay. Each outburst has been 
assigned a letter for reference. 



Relative fluence 

Figure 13. Disc radius versus outburst fluence for the seven fitted 
outbursts of 4U 1705-44. 



and Fig. 13 shows a positive correlation between outburst 
fluence and Rdi B c(L t ). 

If at least some of the knees marked in Fig. 12 actually 
correspond to the disc radius mechanism being considered 
(i.e. they are brinks), then substantial variation in the radius 
is seen. Furthermore, the limit of the exponential decay is 
clearly in several cases below the mean long term count rate 
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Figure 14. The relation between the radius as determined from 
the brink luminosity against that based on the exponential time- 
scale. Circles mark neutron star systems, and triangles black 
holes; individual points can be identified using Table 3. Crosses 
mark decays of 4U 1705-44 other than a and d, as described in 
the text. A diagonal, dotted line marks Rdisc(Lt) = -Rdisc( r e). 

of 12 counts per second, based on the numerical integration 
of the ASM light curve. We account for this by the additional 
theory given in section 5. 

4 SOURCE DATA AND BEST-FITTING 
PARAMETERS 

To force the disc radii based on r e to lie within the range 
ficirc < i?disc(T"e) < b for the two best described systems, 
SAX J 1808.4-3658 and XTE J 0929-314, the viscosity of 
j^kr = 4 x 10 10 m 2 s _1 was adopted, close to the value of 
10 11 m 2 s _1 proposed by KR. To make Rdisc(Lt) consistent 
with -Rdisc(Te) for most of the neutron star systems a value of 
$ = 1.3 x 10~ 12 m 2 s J -1 was adopted. A better fit was found 
for the black hole systems by reducing $ by a factor of ten, 
consistent with the extra factor of Ho discussed in the intro- 
duction. $ is somewhat small compared with the range of 
4 < <E>/(10~ 12 m 2 sJ _1 ) < 9 adopted by KR. Furthermore, 
since -Rdisc is based on $L t , our systematic underestimate 
of L t by considering only the flux in the 1.5-12 keV band 
suggests a yet smaller value of Calculated disc radii using 
these parameter values are given in Table 3. 

Fig. 14 shows the relation between the values of i?disc 
calculated from L t and r e . We note that for 4U 1705-44 the 
results would be more consistent if we adopted the lower 
value of 3>, appropriate to black hole systems. This may 
reflect differences in the source spectrum as well as intrinsic 
differences in $ since we have assumed the same relation 
between 1.5-12 keV count rate and total X-ray flux for all 
systems. In addition to this effect, we note that the results 
for 4U 1705-44 in Fig. 14 do not lie on a straight line. 

To examine this issue we first determined a best-fitting 
value of 3>/^kr from decays a and d, which we deem the 
strongest results by inspection of the light curve. We then 
demanded that this value of $/ikr apply to each decay of 
4U 1705-44. In Fig. 15 we show the effect this has on the fits, 
plotting the residuals for the individual best-fits, and for the 
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Figure 15. Residuals to the fitted portions of the RXTE ASM 
light curve of 4U 1705-44. Left column: using best-fitting pa- 
rameters for each individual decay. Right column: Relating the 
brink luminosity to r e by Nt/r e = 2.8 X 10~ 5 ASM counts s~ 2 
in all cases. 



constrained fits. If the original fits were strongly determined 
by the data, we would expect the residuals to appear signifi- 
cantly worse in the constrained fits. In contrast, we find that 
in some cases the residuals appear to improve, and in none 
is the fit substantially worse; the reduced \ 2 f° r trie whole 
set of fits increases from 11.4 when L t and r e are optimized 
individually to 12.1 for the constrained fitting. We note that 
technique is only applicable when comparing decays of the 
same source, where variation in the quantity &/iskr is min- 
imized. 

The large value of R disc (Lt) for XTE J 1751-305 may 
be again the result of spectral differences, or in this case 
may be the result of an erroneous source distance estimate. 
XTE J 1751-305 has been assumed to lie at a distance of 
8kpc, being possibly associated with the Galactic Centre 
(e.g. Gierlihski & Poutanen 2005); if the true distance is 
more like 3 kpc then the two values of -Rdisc will be approx- 
imately equal. 

The increasing value of -Rdisc(T e ) for the 2002 decay of 
SAX J 1808.4-3658 and the erratically increasing luminosity 
between days 65 and 70 of Fig. 3 are both consistent with 
the interpretation that the outer disc edge was not in the 
hot state at the beginning of the decay. As such, the disc 
radius measured corresponds at first to a hot disc smaller 
than the disc as a whole, which is therefore depleted on a 
shorter timescale. At some point the outer disc is heated and 
enters the hot state. As well as increasing the timescale of 
the decay, this introduces extra material into the hot disc, 
causing rebrightening. 
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Source Porbit [min] ^spin [ ms ] a [10® m] A/2 [-^©] Distance [kpc] 



XTE J 0929-314 


43.6 a 


5.4 a 


3.2 


0.008-0.03" 


5 a 


NS 


4U 1543-475 


1617 ± 12 6 




73±4 C 


2.7±1.0 C 


7.5±1.0 C 


BH 


XTE J 1550-564 


2235 ± 14 d 




86-92 d 


1.4±0.5 d 


2.8-7.6 d 


BH 


GRO J 1655-40 


3775.1 ±0.2 e 




117 ± l e 


2.34±0.12 e 


3.2 e 


BH 


4U 1705-44 










7.3 ± 1.0* 


NS 


GRO J 1744-28 


170009 




1909 


0.2-0.79 


8.5 h 


NS 


XTE J 1751-305 


42^ 


2.3' 


3.1 


0.013-0.035* 




NS 


XTE J 1807-294 


40 J 


5.25 j 


3.0 


0.01-0.022 fe 




NS 


SAX J 1808.4-3658 


120.8 ! 


2.5 1 


6.3 


0.04-0.1' 


2.5 m 


NS 


GX 339-4 


2500™ 




80-120" 


< 1.1" 


15" 


BH 



a: Galloway et al. (2002) f: Haberl & Titarchuk (1995) k: Falanga et al. (2005) 

b: Orosz et al. (1998) g: Rappaport & Joss (1997) 1: Chakrabarty & Morgan (1998) 

c: Park ct al. (2004) h: Nishiuchi ot al. (1999) m: in 't Zand ot al. (2001) 

d: Orosz ot al. (2002) i: Markwardt ot al. (2002) n: Hyncs ct al. (2004) 

o: Orosz & Bailyn (1997) j: Campana et al. (2003) 



Table 1. System parameters from the literature. Orbital separation a is based on a 1.4Mq neutron star where an accretor mass is not 
available. 



Decay L t [ 10 27 Js" 1 ] L e [ 10 27 Js" 1 ] r e [ days ] n [ days ] 



XTE J 0929-314 


49 ±1 


48 ±1 


6.9 ±0.4 


19.7 ± 1.2 


4U 1543-475 


27100 ± 900 




12.05 ±0.13 




XTE J 1550-564 


11200 ± 100 








GRO J 1655-40 


2900 ± 20 


2610 ±50 


13.7 ±0.9 


19.3 ±0.3 


4U 1705-44 a 


2400 ± 300 


1920 ± 80 


25.1 ± 1.2 


22.7 ± 1.4 


4U 1705-44 d 


1680 ± 50 


900 ± 200 


16.9 ± 1.9 


10.0 ±0.6 


4U 1705-44 e 


1650 ± 80 


0-220 


36 ±3 


8.4 ±0.8 


4U 1705-44 h 


3940 ± 40 


3900 ± 50 


12.0 ±1.4 


19.8 ± 1.7 


4U 1705-44 i 


4910 ±20 


4910 ± 20 


19.7 ± 1.0 


36.4 ± 1.6 


4U 1705-44 j 


1620 ± 40 


4170 ± 30 


26.8 ± 1.4 


12.0 ± 1.7 


4U 1705-44 k 


3430 ± 20 


3410 ± 20 


15.9 ±0.7 


50.0 ±0.9 


GRO J 1744-28 1996 


8050 ± 100 


0-40 


46.61 ±0.12 


21.4 ± 1.0 


GRO J 1744-28 1997 


6120 ± 80 


0-170 


42.9 ±0.8 


18.7 ± 1.2 


XTE J 1751-305 


361 ± 1 


110 ±30 


5.9 ±0.5 


1.69 ±0.02 


SAX J 1808.4-3658 1998 


38.54 ±0.03 


32.3 ±0.2 


4.89 ± 0.06 


2.978 ± 0.007 


SAX J 1808.4-3658 2002a 




50 ±7 


2.8 ±0.3 




SAX J 1808.4-3658 2002b 


45.1 ±0.1 


29.6 ±0.6 


5.0 ±0.1 


3.22 ±0.03 


GX 339-4 


1916 ±90 


900 ± 200 


57 ±3 


27 ±4 



Table 2. Decay light curve parameters. Luminosities are based on the photon index and distance from Table 1 and values of Nfj from 
the same sources, and are stated for a spectral range of 1.5-12 kcV. For consistency with Figs. 7 and 5 we take 1 ASM count s^ 1 to 
equal 31 PCA counts s _1 for our purposes, suggesting that T ~ 2.1. 



5 MODIFIED EXPONENTIAL LIMIT 

We have assumed in the initial model that the limit of the 
exponential decay is equal to the assumed constant mass 
transfer rate from the donor star. This requires that the hot 
disc is capable of supporting the additional mass flux — M2 
at all times and radii. However, the exponential nature of 
the decay itself implies that as the hot disc is depleted in 
surface density, its mass transfer rate falls. The mass trans- 
fer equation of the hot disc is linear in surface density, which 
can therefore be decomposed into two components; one rep- 
resenting the decay analogous to an unfed disc, the other 
corresponding to the effects of feeding the hot disc at the 
rate — M2. The former component we denote S e and the lat- 
ter E2. Because the viscosity of the accretion disc depends 
on temperature, and the temperature is declining during the 
decay, it is necessary for £2 to increase with time if it is to 



sustain the same mass flux, implying that the corresponding 
flux fj,2 must decrease inwards. Thus the limit of the central 
mass accretion rate does not correspond to —M2, but rather 
/i2(i?in). The relevance of this issue is made clear by the light 
curves of 4U 1705-44, in which the limit of the exponential 
decay is often below the long term mean count rate of 12 
ASM counts s _1 . 

The standard derivation of the diffusion equation (e.g. 
Pringle 1981) can be used to give the mass transfer rate 
inward through the disc /i as 

V^QkR 1 ' 2 ^^ 2 )' . (25) 
Using the relation 

At' = 2ttR±, (26) 
we derive equation 9 of King (1998) as expected. We require 
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Decay 


-f^circ 


bi 


Rdisc(Lt) 


-Rdisc ( T e) 


Relative 


-Mi 




[ 10 8 m ] 


[ 10 8 m ] 


[ 10 s m ] 


[ 10 8 III ] 


fluence 


[ 10- 12 Af Q yr- 1 ] 


XTF, 10929-314 


1.5—1.9 


2 63—2 82 


o m + 03 


2 fi7 + 08 

i.U 1 _1_ W.UO 




84 ± 2 


4TT 1 ^43-475 


12.7-17.1 


43.4—48.7 


iss + n 5 

1U.O _l_ KJ • <J 


*? 1 fi + 02 

O. _1_ U.Uii 






XTF, 1 1 ^50-564 


20 0-2fi 4 


58 5—64 2 


1 2 1 5 + 05 

lii. lu _l_ \) .KJiJ 








ORO 1 1 65^-40 


21 S— 22 1 


70.9—72.1 


a i A -u n (12 


3.43 ± 0.11 




Afinn + 1 nn 

'-tVJVJVJ _l_ 1UU 


4TT 1 70^ 44 a 

1U 1 I U J 11 a 






17.7 ± 1.1 


ij. J-U U. J- Zr 


0.57 


qaoo 4- 9fin 


4TT 1 70^ AA h 










0.32 




4TT 1 70^ AA r 










1.00 




ATT i 70^-44 H 






14.8 ± 0.2 


4.2 ± 0.2 


0.34 


1UUU _l_ t:UW 


ATT 1 70^-44 p 






14.6 ± 0.4 


6.1 ± 0.3 


0.21 


< 400 


ATT 1 7(1^-44 f 










0.21 




4TT 1 705-44 p- 










0.30 




ATT 1 70^-44 ti 

rr I.J 1 1 \J 1 1 11 






22 6 + 1 2 

jLjZj .\J _|_ U. liJ 


3.5 ± 0.2 


0.69 


kqoo + 200 


4TT 1 70^ 44 i 

1VJ X / UU 11 1 






ZO.O m u.uo 


A cro _|_ n 1 9 

1. HI U. 1 ^ 


n 7fi 




4U 1705-44 j 






14.5 ± 0.18 


5.27 ±0.14 


0.33 


7300 ± 50 


4U 1705-44 k 






21.1 ± 0.06 


4.06 ± 0.09 


0.46 


6000 ± 50 


4U 1705-44 1 










0.29 




GRO J 1744-28 1996 


32.5-52.3 


112.1-133.9 


32.3 ±0.2 


6.95 ±0.01 




£ 70 


GRO J 1744-28 1997 


32.5-52.3 


112.1-133.9 


28.2 ±0.2 


6.67 ±0.06 




< 300 


XTE J 1751-305 


1.4-1.7 


2.5-2.7 


6.85 ± 0.01 


2.47 ±0.11 




200 ± 50 


SAX J 1808.4-3658 1998 


2.06-2.71 


4.68-5.07 


2.238 ± 0.001 


2.252 ±0.014 




57 ±4 


SAX J 1808.4-3658 2002a 


2.06-2.71 


4.68-5.07 




1.70 ±0.09 




90 ± 10 


SAX J 1808.4-3658 2002b 


2.06-2.71 


4.68-5.07 


2.421 ± 0.003 


2.28 ±0.02 




51.8 ± 1.1 


GX 339-4 


> 23 


> 66 


4.99 ±0.12 


6.88 ±0.23 




1600 ± 400 



Table 3. System radii and mass transfer rate. The two values of -Rdisc an d —M2 are calculated from the X-ray light curves; -R c i rc and 
b\ are calculated from the data in Table 1. Errors exclude the effect of the uncertainty in source distance. 



the form given in (25) to prevent a constant of integration 
arising. 

To make use of (25) we must adopt an appropriate form 
for the viscosity v. We use the a- viscosity form 



v — ac s H, 



(27) 



where c s is the sound speed, and the scale height is given by 
equation 2 of King (1998); 



1/2 



(28) 



Using (1) to substitute for T and assuming that H(R) is a 
power law we obtain 



H 



/ k fW (l-A)(n-l) ^ 7 
\GNhrn) \ 4tt<7 



4 /7 r 9/7 , 



(29) 



for a mean particle mass m ~ m p . We write (29) in terms 
of constant Ho as 



H = H(,R disc ("^"^ 



f(£) 



(30) 



i / \ -Rdisc I 

where n = 9/7 is derived and Ho ~ 0.2 is adopted from KR. 
Burderi, King & Szuszkiewicz (1998) note that if the X-ray 
flux is assumed to be absorbed at the mid-plane, the power 
law index is changed to 45/38; we therefore consider values of 
45/38 < n ^ 9/7 as valid to simplify subsequent derivation. 
Using this form for H in (1) to derive the viscosity we find 
that 



(Lx\ 2/7 ( r y 9 "- 3 )/ 8 

\L t ) \R disc ) 



V = Vo 

where uq is defined as 



(31) 



v = a 



k A L t (1 - A) (n - 1) 
4nam 4 



1/8 



9/8 r>3/4 



N H d 



(32) 



We simplify (31) by adopting n = 11/9, so (25) becomes 
6nv ( L x \ 2/7 1/2 



M = 



If we impose fi 
-Mi (L 



"i-KVo 



-M2 for all radii we obtain 



(33) 



(34) 



corresponding to a total contribution to mass in the hot disc 
of 



m2 



-2M 2 
3fo 



-2/7 



(35) 



where the radius of the inner disc edge is much smaller than 
Rdisc- As Z/x decreases, a greater surface density and total 
mass is required to transport the same — M2 inward, so some 
of this material must be absorbed to increase the surface 
density. For a typical luminosity at the beginning of the 
exponential phase 2.5 times greater than L t the change in 
mass required is 

(l-2.5- 2/7 ). (36) 



Am,2 



M 2 



3^o 



Examining 4U 1705-44, for which we have the largest mean 
luminosity relative to L t , we find a mean long term count 
rate of 12 ASM counts s _1 , whereas the limit of the decay 
a for instance is 8.2 ASM counts s _1 . Taken over an expo- 
nential phase duration of 57 days this makes the l.h.s. of 
(36) equivalent to 1.9 x 10 7 ASM units of mass. Using the 
conversion factor from WebPIMMS of 1 ASM count s _1 = 
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4 x 10~ 13 Js _1 m~ 2 , at a distance of 7.3 kpc this integrated 
flux is equivalent to 3x 10 20 kg. Using f KR = 4x 10 1() mT 2 s _1 
(as in Table 3) we find 

R disc ~ 6.4 x fO 8 m, (37) 

in strong agreement with the values given in Table 3 despite 
the weakness of the derivation. In principle it is necessary 
to simultaneously derive the time and radius dependencies 
of fj, given the outer boundary condition of fi(Rdiec) = ~M 2 , 
but this is beyond the scope of this paper. 
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6 CONCLUSIONS 

A knee in the light curve of the decay from outburst of an 
transient LMXB is a natural consequence of mass transfer 
onto the outer edge of the disc, since this supply is effec- 
tively cut off from the compact object when the outer disc 
enters the cool low-viscosity state. When the knee can be 
interpreted in this way we refer to it as a brink, since knees 
of other types may be seen in the decay outburst. The X- 
ray luminosity at which the brink occurs is that at which 
the outer disc edge is just kept hot by central illumination, 
allowing this radius to be calculated. In addition, the ex- 
ponential time-scale of the decay gives a second measure of 
the disc radius; these two estimates are in good agreement. 
By examining systems with well constrained disc radii we 
deduced values of the constants $ and vkr required to per- 
form the calculations. This allows more accurate results to 
be obtained in other systems. Where the source distance is 
poorly constrained, the timescale T e allows the absolute lu- 
minosity L t to be estimated, providing a measure of source 
distance. 

As well as 'secondary maxima' associated with the 
brink, we have identified maxima whose natural explana- 
tion is the discontinuous increase in the radius of the hot 
disc during outburst. This increase arises from initially self- 
shadowed disc regions becoming irradiated. This process can 
occur until the whole disc is in the hot state. It is character- 
istic for the time-scale of the exponential decay to increase 
with each such maximum, and it is expected that the limit 
of each exponential decay will also increase. Although we do 
not have sufficient data to examine the matter in detail, it 
is probable that some secondary 'maxima' of this type will 
also involve monotonic decay, and therefore probably will 
appear as a knee in the X-ray light curve. 

Current instrumentation is capable of measuring the 
light curves of X-ray transients in M31, for example Tru- 
dolyubov, Priedhorsky & Cordova (2006) show part of a 
2004 July decay of XMM J 004315.5+412440. While this 
light curve, with only 4 data points over 3 consecutive days, 
is too sparse to allow us to deduce any parameters from 
fitting the decay, our method could be applied to better- 
sampled extragalactic transient decays. Since the distance 
is relatively well-known for objects in external galaxies, we 
could use our method to deduce their accretion disc radii. 
This would contribute to constraints on the fundamental 
system parameters, such as orbital separation and mass ra- 
tio, which might otherwise be impossible to determine. 
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